Two Phase Transitions in Chiral Gross-Neveu Model in 2 + e Dimensions at Low N 
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We show that the chiral Gross-Neveu model in 2 + e dimensions has for a small number N of 
fermions two phase transitions corresponding to pair formation and pair condensation. In the first 
transition, fermions and antifermions acquire spontaneously a mass and are bound to pairs which 
behave like a Bose liquid in a chirally symmetric state. In the second transition, the Bose liquid 
condenses into a coherent state which breaks chiral symmetry. This suggests the possibility that in 
particle physics, the generation of quark masses may also happen separately from the breakdown of 
chiral symmetry. 



I. INTRODUCTION 
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The Nambu-Jona-Lasino (NJL) model M and its A^-component version, the Gross-Neveu (GN) model ||, are field 
Oh theories of zero-mass fermions with quartic interaction which provide us with considerable insight into the mechanisms 
of spontaneous symmetry breakdown. Both models can formally be turned into pure boson theories. In an SU(3)x 
SU(3)-symmetric version, the NJL model has has been shown to be equivalent to a a chirally SU(3)xSU(3) invariant 
\Q a- model which reproduces all well-known relations of current algebra (|] . For recent work and citations see Q . 

The Gross-Neveu model is exactly solvable in the limit of N — > 00. For an attractive sign of the interaction, a 
collective fermion-antifermion field acquires a nonzero vacuum expectation, and the system shows quasi-long-range 
order. In D = 2 + e dimensions, the order of this state becomes proper long-range. The ordered state is reached in a 
second-order phase transition from a disordered state if the renormalized coupling constant g becomes larger that a 
critical value g* = ire. The disordered state at small g < g* consists of massless interacting fermions. It exhibits chiral 
symmetry, in which fermions are transformed by a phase rotation containing a 75 matrix. In the ordered state at 
' larger g > g* , however, the fermions acquire spontaneously a mass, and the chiral symmetry is broken spontaneously. 
The purpose of this note is to point out, that in a modified Gross-Neveu Model in which pairs of fermions form bound 
pair states analogous to the Cooper pairs in superconductor, the paired phase becomes incoherent in a Kosterlitz- 
Thouless-like transition || if the number of field components N drops below a certain critical value N c < 8. From 
this we conclude that in 2 + e dimensions, the transition in which the pairs form exists independently of a transition 
"*7* ' in which they condense. 

Oh, In the ordinary Gross-Neveu model, the role of the Kosterlitz-Thouless-like transition is played by an Ising transition, 
■ which appears in addition to the transition in which the collective state forms. 
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The original Gross-Neveu model has the following 0(AT)-symmetric Lagrange density 

C = $ai$ll>a + {i'a^a) 2 ■ (1) 

where the index a runs from 1 to N . At the mean-field level, the effective action is equal to the initial action 

r *]=a $] = J d D x [m$* + ^ (* Q *a) 2 ] • (2) 

In general, we obtain all Green functions from the generating functional 

Z[rj,ff\ = e* 1 ' 1 ' 1 = / V^ViPe lA ^^ + ' l ^ ri+c - c -\ (3) 



II. PROPER GROSS-NEVEU MODEL 
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where rj(x) and T)(x) are fermionic anticommuting sources. A collective field a ~ gip^P is introduced to rewrite (0) as 



The fields ^>(x) are integrated out according to the rule to yield a generating functional containing only the collective 
field a(x): 

Z[r),fj] = J p^e^oiiM-^^ ( 5 ) 

with the collective action 

AouM = iV |-^-<7 2 -iTr log[^ -tr(a:)]J. (6) 

where Tr denotes the functional trace. 

In the limit N — > oc, the field cr is squeezed into the extremum of the action and we obtain the effective action 

ir [e, ^ = "^ e2 ( x ) ~ iTr lo s W ~ E W1 + [*0 - E 0*0] *« ( 7 ) 

The extremum of T [E, ^1 is given by the equtions of motion, 

[i0 - H(x)] *„(x) = 0, (8) 

1 - i 

E(at) =g- trGs(a:,a;) - — goip a ip a (x), G^(x,y) = — -, (9) 

Jy lip — ia 

where the trace symbol tr is restricted to the Dirac indices. The expectation a of a fermionic field is always zero, 
so that we only must solve the gap equation 

X(x)=g trGj;{x,x). (10) 
Thus, as far as the extremum is concerned, we may study only the purely collective part of the exact action 

l r [E] = --^E 2 - i TrlogzG E 1 . (11) 

^ £go 

The ground state is given by a constant gap field Eo, for which ( |l0|) yields either Eo = or 

where we have performed a Wick rotation p° — > ip A to euclidean momenta p^ = (p 1 , p 2 , p 3 , p A ) with the metric 
p\ = —p 2 . The Dirac matrices have dropped out, except for the unit matrix whose trace is 2 D / 2 for even D. This 
expression may be continued to any non-integer value of D. 

For a constant E, the effective action gives rise to an effective potential 



-«(S) = --T[E] =— E 2 - tr (1)- [ log \p\ + E 2 1 

N V ' N 2go y '2 J (2ir) D B lFE J 
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Performing the integral yields in D = 2 + e dimensions with e > 
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E 2 /E^ 2+£ 
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(13) 



(14) 



where \x is an arbitrary mass scale, and the constant b e stands for 



b c = -V 2 S D lW2)r(l - D/2) = l—L^TQ - D/2), (15) 
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which has an e-expansion b e 
introduced by the equation 



so that 



[1 — (e/2) log (2ne 7 )]/7re + 0(e). A renormalized coupling constant g may be 



be = 



1 



Extremizing this we obtain either So = or a nonzero So solving the gap equation (|12|) in the form 

9* D /S S ' 



1 



/' 



(16) 



(17) 



(18) 



where g* = —l/b e ~ 7re. A nontrivial solution of this is called gap. It specifies the mass which the fermions acquire 
from the attractive interactions, and will be denoted by M. The second derivative of w(S) shows that the solutions 
So = and So ^ are stable for go > and go < 0, respectively. Denoting the solution So of (18) for g = oo by 
Mqo ■ we may write the g-dependence of the gap as 



M{g) = Moo 1 - — 
V 9 

In terms of M, the effective potential (O) can be rewritten as 
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It has a minimum at S = M, where it yields the condensation energy v(M) = —NM D /4:TT. 



(19) 



(20) 



III. CORRELATION FUNCTIONS OF PAIR FIELD 



If N is no longer infinite, the pair field a in the partition function (^) performs fluctuations around the extremal 
value So = M. For large N, the correlation functions of a(x) can be extracted from the leading effective action (Q) 
at Sq = M. Setting S(x) = M + S'(x), we expand 



implying a propagator of the c'-field 
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where n(<7) is given by the self-energy diagram 

U(q) - ~2 D ' 2 



d D k 



N 1/go + U(q) ' 

k(k-q) E -M 2 



(21) 



(22) 



(2n) D (fc| + M 2 )[(fc- g )|+M 2 ]' 



th integral being performed over euclidean D-dimensional energy-momenta. With standard Feynman methods, this 
can be transformed into a simple integral 



n(«) 



D(D-1) 



b f M e I dx 
o 



(23) 



Inserting here Eq. (18), we obtain 
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- + n( 9 ) = // [\-- 

.90 \g* g 



which can be expanded for small q as 

1 



-1 + (D-1) / dx 
Jo 



M 2 



x(l - x) + 1 



e/2' 



-+U( q )=e,f (1-1 



l+P-i; 
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12 M 2 



This shows that the propagator ( |22| ) has a correlation length 



Inserting the g-dependence of M from (p_9|), we see that 



Moo V 12 / V 3 
so that the coherence length diverges for g — > g* with a critical exponent f = 1/e. 



(24) 



(25) 



(26) 
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A. Chiral and Complex Pair Field Version of Model with Goldstone Bosons 



We want to prove the existence of two phase transitions in the chiral version of the Gross-Neveu model, whose 
Lagrange density is 



go 

2N 



The collective field action (^J) is then replaced by 

N 



AoiiM = N \ --— {a 2 + it 2 ) - iTr log [»0 - a{x) - i 75 7r] 

%go 



(28) 



(29) 



This model is invariant under the continuous set of chiral 0(2) transformations which rotate a and it fields into each 
other. This model is equivalent to yet another one which is closely related to the BCS model of superconductivity. 
Its Lagrange density is 



Here C is the matrix of charge conjugation which is defined by 



C'fC 



(30) 



(31) 



In two dimensions, we choose the 7-matrices as 7 = a 1 , 7 1 = —ia 2 , and C = 7 1 . Note that (^ a C^J)^ = %pJCtp a , 
implying that go < corresponds to an attractive potential. The second model goes over into the first by replacing 
ip — > |(1 — 75)^ + 4(1 + 75)2C-0 T , where superscript T denotes transposition. In the Lagrange density (30) we 



introduce a complex collective field by adding a term (N/2go) A — jj-ipf;Cipb\ , leading to the partition function 



Z[r],fj] = / VipVipVAexp 



$ a iHa + \ (Aty^Cty* + c.c.) + iPr) + # - ^|A| 5 
2 2go 



(32) 



The relation with the previous collective fields a and n is A = a + Z7r. In order to integrate out the Fermi fields we 
rewrite the free part of Lagrange density in the matrix form 



)(<& 



(33) 
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which is the same as ipifitp, since ip T CCip T — tptp, ip T C <J) Cip T = iptflip. But then the interaction with A can be 
combined with ( |33| ) in the form ^tpfG^cf), where 



_ ( i> 



, </? = (f,4ic- 1 ) (34) 



are doubled fermion fields, and 



«X-{ C oc) (*%)--<&)' <*> 

is the inverse propagator in the presence of the external field A. Now we perform the functional integral over the 
fermion fields, and obtain 

Z\j] = J VAVA^e mA ^+^ G ^ , (36) 



where A\A] is the collective action 



^[A] = -i|A| 2 -^Trlog J G A 1 (37) 



and j a is the doubled version of the external source 

(38) 

This is chosen so that $r] + fjip = \ {j T 4> — 4> T j)- 111 the- limit N — > oo, we obtain from ( |36"| ) the effective action 

—T\A, *1 = — |A| 2 - -Tr logiGT 1 + - 
N 2go 2 N 



-r[A, *] =— |A| 2 - -Tr logzG^ 1 + -tf.iGk 1 *, (39) 



in the same way as in the last chapter for the simpler model with a real a-field. 

The ground state has \& = 0, so that the minimum of the effective action implies for Ao either Ao = or the gap 
equation 

l = fTrG Ao , (40) 
where we may assume Ao to be real. With the Green function 

the gap equation (^0|) takes the same form as (|l2|): 

/ s f d D p 1 

where we have again set M = A$. The renormalization of the coupling constant and of the effective potential yields 
the same equation for v(Aq) — v(M) as before, so that the previous stability discussion for g < g* and g > g* holds 
also here. 

Let us now study the propagator of the complete A-field. For small deviations A' = A — Ao away from the ground 
state value we find from ( |39| ) the quadratic term 



1 -6 2 T= --4 — -Tr 



A't \ „ ( A't 



N 2 t 5o 2 

The second term in curly brackets may be written more explicitly as 



A' j Gm [ A' I Gm 
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and becomes 



\ {A/ 2 (A' 2 + A't 2 ) U(q 2 E /M 2 ) + 2|A'| 2 [l%|/Af 2 ) - M 2 fl(q 2 E /M 2 )] } , 



where H(q% / M 2 ) is the previous self-energy (|23|), and II (q%/M 2 ) is the function 



ft (q 2 E /M 2 ) = i2 D ' 2 



d k i i D, ,„ ,„ N f , n r 9 .„ . , r9n Z5/2-2 



As a result, the action for the quadratic deviations from Aq = M can be written as 



-+A ] |A 

5o 



,'|2 



^B(A 



with the coefficients 



12 | A'*2\ 



A = --b t M e [{D - l)J{ (q%/M 2 ) - (D/2 - l)JI(q 2 E /M 2 )] , B = -b e M<(D/2 - 1) J 2 e (q%/M 



and the integrals 



d D x [zx{l - x) + l] D ' 2 ~ x , J e 2 {z)= d D x[zx(l-x) + l] D/2 - 2 



Jt{z) - 

Thus the propagators of real and imaginary parts of the field A' are 

i 1 „ i 



Ga 



> Ia A' TC 



Ga< 
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Nl/go + A + B' ~ A ™ A ™ Nl/go + A-B' 
The excitation spectrum is given by the zeros of the denominator functions 

— + A + B = ^b e M e [1 - (X> - l)Ji (q%/M 2 )] , 

-+A-B= ^b e M e { [1 - (D - 1)J{ (q 2 E /M 2 )] +{D- 2) J\ (q 2 E /M 2 ) } 

go z ~ 

By expanding J[(z), J^iz) in powers of z = q E /M 2 « 0, 

•Kit) ~ 1 + + 0(z 2 ), Jl(z) ~ 1 + ^—^z + 0(z 2 ), 



we find 
1 



D, 



9o 



12 



D - 1 



12 



A + B = b e M e (D - 2) 1 + z + 0{z 2 ), —+A-B 



12 



,9o 



2 e 12 



3z + 0{z 2 



Inserting here the gap equation (JlOj) , we obtain 



50 



-A+B = 



1 1 

g* g 



1 + R=lA 

12 M 2 
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-A-B = 



g* gj 4M 2 



(43) 

(44) 

(45) 

(46) 

(47) 

(48) 
(49) 

(50) 
(51) 

(52) 



Recalling (g5[) we see that the propagator of AJ. C coincides with that of £' in the standard Gross-Neveu model, so that 
the fluctuations of AJ. C have the same correlation length (26), with at a critical exponent v in ( p7| ) as g approaches g* . 
In contrast, the propagator of the imaginary part of ATas now a pole at q 2 = 0: 



14/1 1 

<-■&'.&'. — T7- — 

un N e \g* g 



M — + regular part at q = 0. 



(53) 



The sign of the pole term guarantees a positive norm of the corresponding particle state in the Hilbert space. The 
particle is a Nambu-Goldstone boson. 
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IV. SECOND PHASE TRANSITION 



We are now prepared to show that the pair version of the chiral Gross-Neveu model in 2 + e dimensions has two 
phase transitions. Consider first the case e = where the collective field theory consists of complex field A with 
0(2)-symmetry A — > e^A. From the work of Kosterlitz and Thouless || we know that such a field system possesses 
macroscopic excitations of the form of vortices and antivortices. These attract each other by a logarithmic Coulomb 
potential, just like a gas of electrons an positrons in two dimensions. At low temperatures, the vortices and antivortices 
form bound pairs. The grand-canonical ensemble of pairs exhibits quasi- long-range correlations. At some temperature 
T c , the vortex pairs break up, and the correlations becomes short-range. The phase transition is of infinite order. 

This transition is most easily understood in a model field theory involving of a pure phase field 9(x), with a Lagrange 
density 

C = ^[d9(x)]\ (54) 

where j3 is the stiffness of the ^-fluctuations. The important feature of the phase field 9 is that it is a cyclic field with 
9 = 9 + 27r. In order to ensure that such jumps by 2ir carry no energy, the gradient in the Lagrange density needs a 
modification which allows for the existence of vortices and antivortices. This will not be discussed here in detail, since 
the reader may consult the literature for it [HQ] . We only state here that after including vortices and antivortices at 
positions x%, Xj, their partition function can be written as 



^2 exp i 4tt 2 /3 Uli lo S 

gas ^ i<j 



(\xi - Xj\/r ) | , (55) 

where rg is the size of the vortices. For a single vortex-antivortex pair, the average square distance r 2 is 

< r 2 >cx / drrr 2 e~ 27r/31og ( r / ro ) oc . (56) 

This diverges as the stiffness falls below /3kt = 2/ir k 0.63662. A more detailed study shows that this is an exact 
result for a very dilute system of vortices and antivortices . 

The large-stiffness state with bound vortex pairs has a coherent phase field 9(x), the low-stiffness state with 
separated vortex pairs exhibits incoherent phase fluctuations. The same situation is found in three dimensions, only 
that the excitations are vortex lines. These become infinitely long and prolific in a second-order phase transition Q] 
at a critical point [3 C ~ 0.33. 

The result ([33|) for e = can now be used to estimate a critical value of the number of field components N = N c 
below which the phase fluctuations of the complex field A' become so violent that the system has a phase transition. 
For this we write AJ m = M9 and find from ( |53| ) a propagator of the 0-field 

1 4:7T 

Gee ~ T7^7 + regular terms. (57) 
N q z 

Comparing this with the propagator for the model Lagrange density (|5 

1 i 

PV 2 



Gee = -- (58) 



we identify the stiffness (3 = N/Att. The pair version of the chiral Gross-Neveu model has therefore a vortex-antivortex 
pair breaking transition if N falls below the critical value N c = 8. 

Consider now the model in 2 + e dimensions where pairs form at g = g* rj 7re. A comparison between the propagator 
(|53|) and fl5S| ) yields a stiffness of phase fluctuations 



Air \ g 

The linear vanishing of the stiffness with the distance of the coupling constant g from the critical value g* is in 
agreement with a general scaling relation, according to which the critical exponent of bending rigidity should be equal 
to (D — 2)v. Since the model has v = 1/e [see (p7|)], this yields (D — 2)v = 1, which is precisely the exponent in 
Eq. ©. 
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FIG. 1. The two transition lines in the N — (/-plane of the chiral Gross-Neveu model in 2 + e dimensions. For e = 0, the 
vertical transition line coincides with the iV-axis, and the solid hyperbola degenerates into a horizontal line at N c = 8. The 
quark masses and chiral properties are indicated. 



The stiffness ( |59| ) implies the existence of a phase transition in the neighborhood of two and in three dimensions at 
roughly 

N c w 8 ^1 - tjj , fl«2, N c w 4.19 (l - , D = 3. (60) 

As AT is lowered below these critical values, the phase fluctuations of the pair held A become incoherent and the 
pair condensate dissolves. The different phases are indicated in Fig. El In the chral formulation of the same model, 
the intermediate phase has chiral symmetry in spite of a nonzero spontaneously generated "quark mass" M ^ 0. 
The reason why this is possible is that the "quark mass" depends only on |An.|, thus allowing for arbitrary phase 
fluctuations preserving chiral symmetry. 

The sceptical reader may wonder whether the solid hyperbola in Fig. [I] is not simply the proper (albeit approximate) 
continuation of the vertical line for smaller N. There are two simple counterarguments. One is formal: For infinitesimal 
e the first transition lies precisely at g = g* = 7re for all N , so that the horizontal transition line is clearly distinguished 
from it. The other argument is physical. If N is lowered at some very large g, the binding energy of the pairs 
increases with 1/N {in two dimensions, the binding energy is 4Msin 2 [7r/2(7V — 1)]}. It is then impossible that the 
phase fluctuations on the horizontal branch of the transition line, which are low-energy excitations, unbind the strongly 
bound pairs. This will only happen in the limit N — > oo where the binding energy becomes zero and the two transition 
curves merge into a single curve. This is the situation in the BCS theory of superconductivity, where Cooper pair 
binding and pair condensation coincide. 

In the ordinary Gross-Neveu model, the analog of the phase disordering transition is an Ising transition, in which 
the vacuum expectation value of a jumps between —So and So in a disorderly fashion. In two dimensions, this occurs 
at some critical value N c . In 2 + e dimensions, this transition should again exist independently of the transition at 
which the system enters into a state of nonzero So. It will be interesting to see these two transitions in either model 
confirmed by Monte-Carlo simulations. 
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